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Abstract: We present a model of two-kinks resulting from an explicit composition of two 
standards kinks of the model based on the procedure of Ref. [1]. The two-kinks have 
an additional parameter accounting for the separation of the standard kinks of model. 
We have shown that the two-kinks have two discrete internal modes besides the zeroth 
mode and the continuous spectrum. This new feature signalizes that the head-on collision 
a two-kinks/two-antikinks pair exhibits a rich and complex behavior due to the additional 
channel from which the energy of the system can be stored. We have exhibited the fractal 
structure associated with the main configurations after the collision. We have inferred the 
fractality as the imprint of the nonlinear exchange of energy into the two discrete internal 
modes. 
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1 Introduction 

In a very interesting paper, Uchiyama [1] had proposed an extended hadron model based 
on a modification of the sine-Gordon equation. The objective was to explain the quark 
confinement mechanism using two exact solutions: the hrst is for the modified sine-Gordon 
equation of the neutral scalar field, and the second is for the quark field equation in which 
the quark held interacts with the scalar held [1]. The scalar held solution is obtained by an 
explicit composition of two kinks of the sine-Gordon model and for this reason Uchiyama 
named the new conhguration as the two-kinks. 

Some years ago Bazeia et al. [2] introduced a new class of topological defects in systems 
described by real scalar helds in (T*, 1) spacetime dimensions. By restricting to D = 1 they 
have considered the family model with the potential, 

= ( 1 . 1 ) 

where the parameter p is related to the way the scalar held self-inter acts. This model can 
be understood as a generalization of the model, which is recovered for p = 1. For p even, 
the case p = 2 is special and describes an unstable lumplike conhguration [2]. For p = 4,6,.. 
the potential (1.1) requires that </> > 0 or (^ < 0 under the change (j) —<P- As pointed out 

by Refs. [2, 3] it is possible to construct topological defect in the form = tanh^(x/p) 
(x > 0) and = — tanh^(x/p) {x < 0). Moreover, non-topological lumplike defects can 
be envisaged according to Ref. [4]. It can be shown that all defects with p even, topological 
or non-topological, are unstable and therefore of few interest. 

The new static structures arising when p = 3, 5.. connect the minima cp = ±1 passing 
through the symmetric minimum at <p = 0 are called two-kinks defects [2]. As shown in 
Fig. 1, the two-kinks defects seem to be composed of two standard kinks symmetrically 
separated by a distance proportional to p. The energy density prohle reinforces this view. 

The static two-kinks solution associated to the potential (1.1) is described by (p{x) = 
±tanh^(x/p). By implementing standard perturbation about the two- kinks, (p{x,t) = 
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4>ix) + e^‘^^x{x), one can derive straightforwardly the following Schrodinger-like equation 
for the eigenmodes Xnix), 



It happens that the effective potential diverges at the origin also shared by kinks of [5]. In 
other words, the mass squared aX. (j) = Q diverges. This is not good for quantization since it 


cannot be used as a perturbation ground state. 

Two-kinks can be applied to the brane-world scenario [6-8], and in connection with 
processes of their formation in perturbed sine-Gordon models [9]. Another possibility is 


to use two-kinks to describe the magnetic domain walls in constrained geometries [10, 11]. 


Recently, we have studied numerically the head-on collision of two-kinks defects where we 
have noticed the appearance of unstable lump-like structures of non-topological nature as 
critical configurations [12]. 

Although kink structures are ubiquitous in many branches of physics, particular inter¬ 
est relies on representing them as particles. In this context, kink-antikink configurations 
can mimic particle-antiparticle interactions and may provide insights towards fundamental 
concepts like crossing-symmetry in particle physics scattering [13, 14]. 

In this paper, we have extended this construction considering the two-kinks that arises 
from the union of two standard kinks of the (j)'^ model separated by a distance R. In Section 
2, we have obtained the potential corresponding to the new two-kinks and discussed their 
stability. In Section 3 we have solved numerically the Schrodinger-like (1.2) associated to 
the two-kinks to determine the internal modes. We have found that the two-kinks have two 
discrete internal modes besides the zeroth mode. In general, the eigenfrequencies u depend 
on the separation parameter R. We have also shown that the two-kinks have a continuous 
spectrum of modes with eigenfrequencies u > 2 (cf. (2.1)). Motivated by the existence of 
the additional discrete mode, we have discussed the head-on collision of a two-kinks pair 
briefly in Section 3. We have exhibited a fractal structure with respect to the possible 
configurations after the collision. Finally, in Section 4 we summarize the main results and 
present a brief description of the further developments. 

2 Uchyama’s approach and the new two-kinks model 

Uchiyama had considered the sine-Gordon model whose potential is. 



( 2 . 1 ) 


where for the sake of simplicity m = A = 1 so that the sine-Gordon kink solution reads as. 


<^±(3^) = 4 tan ^ e^^^ 


( 2 . 2 ) 


Then, Uchiyama proposed a new kink solution by taking. 


(/>(x) = (/>+(x) - (j)-{x). 
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(2.3) 










Figure 1. Profiles for the kink and its energy density for p = 3 and u = 0.4. The resulting 
structures are composed of two standard kinks symmetrically separated by a distance proportional 
to p [2]. 


The quantity 2R is the distance between the kinks, and R is identified as the radius of the 
bag in the Uchiyama’s model. After a straightforward calculation, it follows that, 


tan 


(p{x) 


tan(0+/4) — tan((;(!)_/4) sinh(x) 


1 + tan((/)+/4) tan((/)_/4) cosh(2i?) 


The potential for the above solution obtained from the Klein-Gordon equation. 


= sin (j) — (cosh R) ^ 
dx^ 


sm( 


' + 2 sin 1 — 


(2.4) 


(2.5) 


We recover the sine-Gordon equation when the two-kinks are far away, i. e., for R —)• oo. 
The potential of the modified sine-Gordon equation becomes. 


l/(^) = 1 — cos (p + (cosh R) ^ 


3 -|- cos 4> + i cos 



( 2 . 6 ) 


The potential is invariant under the symmetry p ^ p and has period dvr. There is a high 
local minimum 8(coshii)“^ at (^ = 0. Uchiyama have chosen R = 5.06770956457570 such 
that the potential has the values 0.0012687 at p = 0 and 2.0031718 at p = tt. Using value 
of R the sine-Gordon and the above potentials almost coincide. 

We can extend the above procedure to construct two-kinks using the standard kink 
solution of the p'^ model. Let us consider the single kink. 


P±{x) = =t tanh(x — xq), 


(2.7) 


where ± indicate the kink/antikink solution. We propose the following two-kink solution: 


p{x) = p+{x) - p-{x) 


sinh(2x) 

cosh(2ii) -k cosh(2x) ’ 


( 2 . 8 ) 


where the parameter R dictates the separation between the single kinks. If i? = 0 we recover 
the standard kink (2.7). From the Klein-Gordon equation, the corresponding potential can 
be calculated straightforwardly resulting in the following expression. 


_ 2(1 - p'^f [l -k (cosh^ 2R - l)p‘^] 

cosh 2R -k Y^l -k (cosh^ 2R — Vjp"^ 


(2.9) 
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It is instructive to expand the potential in power series of i?, 


^ ( 0 " - 1 )" + ( 0 " -lfR^ + 0{R^), ( 2 . 10 ) 


where we recover the potential of the model if i? = 0 as expected. The potential has 
two minima at (;i> = ±1 and depending on the value of R the origin can be a maximum, 
0 < R < R^ = cosh“^(2)/2, or a minimum if ii > R^. For the sake of illustration, we 
present in Fig. 2 the potential for several values of R as well as the scalar field and the 
energy density profiles of the two-kinks. 

We have noticed that the two-kinks (2.10) have already been derived in a distinct 
context. Chumbes and Hott [7] have considered a model with two interacting scalar helds 
in brane world scenario with the introduction of a superpotential. They have constructed 
new models with one scalar field after establishing an orbit equation [8] relating both fields. 
It turned out that one of the models of two-kinks coincides with we have obtained using 
the Uchyiama’s approach. 



Figure 2. (a) Potential curves for R = 0,0.5,1.0 and 3.0, where for R > 0.6584789484624083 the 
origin turns to be a minimum, (b) Profiles of the two-kinks for R = 2.7 and Bazeia’s two-kinks [2] 
with p = 3 (dashed line), (c) Typical energy density profile of a static two-kink with R = 2.5. 


The stability of the two-kinks (2.8) can be inferred using the approach found in Ref. 
[1]. In summary, it consists in perturbing the two-kink configuration by (/> = (j)o{x) + er]{x), 
where (po{x) is the two-kinks (2.10), e <C 1, and r]{x) is an arbitrary differentiable function 
satisfying the boundary conditions r](zizoo) = 0. It can be shown that the total energy. 
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( 2 . 11 ) 



+m = 


— Eq + eJi + J2 + E>{e^), 


where the stability is fulfilled since Ji = 0 and J 2 > 0. 


3 Internal modes of the two-kinks 

We can determine the internal modes of the two-kinks after solving the Schrodinger-like 
equation (1.2). The zeroth mode or ground state cj = 0 has a known analytical solution 
regardless the kink model under consideration. One can show that, 

Xo[x) oc ———, (3.1) 

where (pxix) is given by Eq. (2.10). In order to determine the other internal modes, we 
have to solve numerically Schrodinger-like equation (1.2) to seek for regular solutions. 

At this point we recall that the cj)'^ model {R = 0) has one discrete mode with eigen- 
frequency u! = \/3- This discrete mode is the shape mode that plays a relevant role in 
explaining the bounces found in the kink-antikink collision. There is also a continuous 
spectrum with eigenfrequencies tOk = + 4 known as the boson modes since they be¬ 

come plane waves asymptotically. 

We have integrated the Schrodinger-like equation using a fourth order Runge-Kutta 
algorithm to determine localized and well-behaved solutions. We have set the initial condi¬ 
tion (x(0) = 0,x^(0) = 0.2) and adjust the value of uj to obtain the desired eigensolutions. 
The novelty presented by the two-kinks is the existence of two discrete modes besides the 
zeroth mode. In addition, we have found a continuous spectrum with a; > 2 no matter is 
the value of R. We have summarized the results in Fig. 3 with the values of lo in function of 
some values of R. For the sake of completeness, we have also displayed the results in Table 
1. Notice that by increasing R, the highest discrete eigenfrequency tends to \/3, whereas 
the smallest tends to zero. Therefore, for large R the two-kinks have effectively only one 
discrete mode. We have presented in Fig.4 the plots of the eigenfunctions of the two-kinks 
for R = 2 together with shape mode of the 4>‘^ model. 

We have also considered the two-kinks associated to the potential (1.1), where similarly 
to the standard kink of the (/>^ model, there is only one internal discrete mode. For instance, 
if p = 3 then cj ~ 0.55288762; there is also a continuous spectrum for ui > 2. 


4 A glimpse on the two-kinks pair collision 

Let us present some preliminary result about the collision of a two-kinks pair. The initial 
conhguration is. 


(/)o(0, x) = -1 0x(O, X -h xo) - (pKiO, X - Xo), 


(4.1) 
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Figure 3. Eigenfrequencies Uk corresponding to the internal modes Xk{x) of the two-kink solution 
with respect to R. The dashed region lu > 2 denotes the continuous spectrum, and the dots indicate 
the discrete modes. 




X 


Figure 4. Discrete eigenfunctions of the two-kinks for i? = 2 corresponding to the frequencies 
ui = 0.165313536 (continuous line) and uj = 1.79331783 (dashed line). The blue line accounts for 
the shape mode of the cj)^ where w = ~ 1.73205081. 


R 

iO 

0 

V3 « 1.73205081 

0.5 

1.42413098, « 2.0 

1.0 

0.84029378, 1.97949332 

1.5 

0.39826056, 1.87313903 

2.0 

0.16531354, 1.79331783 

2.5 

0.06398212, 1.75565933 

3.0 

0.02400534, 1.74077421 


Table 1. Single precision values of the discrete eigenfrequencies corresponding to distinct R shown 
in Fig. 3. 

where 
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(l)K{t,x) 


sinh 2^(t, x) 

cosh 2^{t, x) + cosh(2ii) ’ 


(4.2) 


with ^(t,x) = represents the boosted two-kinks with velocity parameter v. We 

have integrated the Klein-Gordon equation with the above initial configuration using the 
spectral code of Ref. [12]. We have set xq = 15.0 to guarantee that the two-kinks are 
sufficiently apart each other before the collision. In the numerical experiments, we have 
monitored the error rms in the energy conservation [12] and verihed that it is not superior 
to 0.01%. In order to achieve this error we have adopted N = 150 in the spectral expansion 
of the scalar field, namely, (j){t,x) = where aj{t) are the unknown modes 

and 'ipj{x) the basis functions (see Ref. [12] for details). 

There are two free parameters, the initial impact velocity and the separation parameter 
R. If we set large i?, the second discrete internal mode is very close to zero, and the two- 
kinks have one discrete mode with w ~ \/3- Nonetheless, the collision turns to be highly 
complex since the system behaves as a collision of four standard kinks of the 4>‘^ model. 
Although we expect an intricate behavior, we will consider the collision of two-kinks with 
large R elsewhere. 

On the other hand, if R is small, say R < 3, both discrete internal modes are effectively 
present. Several works have shown the role played by the shape mode in the collision of kinks 
in the model. Accordingly, this mode can absorb enough energy from the kink-antikink 
system preventing the kinks from moving apart after the collision and several bounces, or 
resonant configurations can occur [15-19]. The bounces are a consequence of a competition 
of energy transfer to the shape and the translational mode. Therefore, the presence of at 
least one vibrational mode is crucial for the bouncing mechanism. 

We have fixed R = 1.5 in all numerical experiments and evolved the system until 
tf = 600. In this case, both discrete modes of the two-kinks might act as channels from 
which energy is transferred. As a consequence, we expect an intricate and complex behavior 
with the presence of multiple bounces and escape windows resulting from the nonlinear 
competition of the energy transfer to the discrete modes. In this vein, we expect some form 
of fractal pattern related to the configurations that emerge after the collision. However, we 
must recall that fractal patterns in the interaction of kinks of the (j)^ model, although these 
kinks do not have a discrete internal mode [20]. 

The numerical experiments have indicated that there exist a critical velocity, Vc ~ 0.29, 
above which the collision is almost elastic, and both two-kinks escape to infinity after a 
brief interaction. Below the critical velocity, the collision is inelastic, therefore losing part 
of the scalar field. We have noticed the formation of a bound oscillating structure at the 
origin mixed with windows of escape. Moreover, multiple bounces may usually occur before 
the formation of an oscillon at the origin or when the two-kinks escape. For the sake of 
illustration, we have set u = 0.23 to show the formation of multiple bounces before a bound 
solution at the origin is settled down. 

We can exhibit the fractal structure related to the main configurations that emerge 
from the two-kinks interaction for v < Vc- These configurations are: (i) the scalar field is 
trapped at the origin forming a bounded oscillating structure; (ii) the two-kinks escape to 
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Figure 5. Three-dimensional plot of the energy density for the collision of a two-kinks pair for 
u = 0.23, R = 1.5 and initial separation xq = 15.0. There are fourteen irregular bounces before a 
bound oscillatory state forms at the origin. 

infinity; and (iii) two-kinks escape to infinity but leaving behind part of the scalar held 
trapped at the origin. In the dynamical system language, we call these states as basins of 
attraction. Then, in order to look up more carefully to the boundaries between these basis of 
attraction we have set the sequence of initial velocities v = 0.5,0.6, ..,0.29, 0.30 and marked 
the corresponding basis of attractions (i), (ii) and (iii), respectively, by black, white and 
gray strips. These strips are centered about each velocity revealing a coarse distribution of 
the basis of attraction (cf. Fig. 6(a)). The next step consists in selecting a smaller interval 
of initial velocity, say from u = 0.13 to u = 0.14 that is indicated by a dashed rectangle 
shown in Fig. 6(a). We have simulated the two-kinks interaction with velocities inside 
this interval, v = 0.131,0.132, ..,0.139, and again placing the corresponding strips (cf. Fig. 
6(b)). We have repeated further this procedure by taking smaller intervals of velocities, 
0.130 > V > 0.131 and 0.1301 > v > 0.1302 and keeping track of the hnal configurations. 
We have presented the results in Figs. 6(c) and 6(d). The successive rehnements establish 
quite clearly the non-smooth or fractal boundaries that separate the basins of attraction. 

5 Further perspectives 

In this paper, we have introduced a new model of two-kinks by extending the Uchyiama 
approach to the standard kinks of (j)^ model. There is a parameter R asserting the separation 
of the two standard kinks that forms the new structure. In this model, the effective potential 
of the Schrodinger-like equation for the internal modes is well behaved, contrary to the case 
of two-kinks with potential (1.1). We have also briefly discussed the stability of the new 
two-kinks. 
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Figure 6. Fractal structure related to the main configurations after the two-kinks interaction for 
V < Vc- The colors black, white and gray denote, respectively the following configurations: formation 
of an oscillon at the origin, escape of two-kinks, and escape of two-kinks together with an oscillon 
at the origin. Figure (a) provides a coarse distribution of these outcomes represented by strips 
centered at the velocities v = 0.05, 0.06,0.38,0.29,0.31. Figure (b) shows the structure within 
the interval 0.130 > u > 0.14 indicated in (a) by a dashed rectangle. The successive refinements 
are present in Figs, (c), (d) showing that the boundaries between the main final configurations are 
fractal. 

We have solved numerically the Schrodinger-like equation (1.2) to determine the in¬ 
ternal modes of the two-kinks. The novelty is the presence of two discrete modes besides 
the zeroth mode (w = 0). We have also exhibited the influence of the parameter R on the 
eigenfrequencies associated to the internal modes. The existence of a continuous spectrum 
for w > 2 is a common feature no matter value of the parameter R. We have also inspected 
numerically the Schrodinger-like equation for the two-kinks of the potential (1.1) [2], and 
we have found just one discrete internal. 

We have studied the collision of two-kinks for R = 1.5, where both discrete modes are 
effectively present. We have found that \i v > Vc with Vc ~ 0.29, the collision is almost 
elastic and the two-kinks escape to infinity. On the other hand, for v < Vc we have exhibited 
a fractal structure associated with the main configurations after the collision. Moreover, 
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multiple bounces can intermediate any of the final configurations (cf. Fig.5). In general, 
possibly due to competition of energy transfer from the system to both internal modes, we 
have observed multiple bounces and windows of escape. 

This first study does not exhaust the complete understanding of the two-kinks collision. 
There are other models of two-kinks whose interaction could produce identical features we 
have observed such as several windows of escape. Further, we can gain analytical insight 
by using the collective coordinates approach in analogy to the studies of the model. In 
particular, the role played by the two discrete modes in the dynamics of the interaction of 
the two-kinks is a subject worth of investigation. 
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